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Abstract 

We compute the 0{as) QCD corrections to the inclusive process e+e” —> tih. Al¬ 
though the total rate is small, it has a distinctive experimental signature and can po¬ 
tentially be used to measure the top quark-Higgs boson Yukawa coupling. The QCD 
corrections increase the rate by a factor of roughly 1.5 for e'’'e“ —> tih at y^=500 GeV 
and Mfi = 100 GeV. At ^/s = 1 TeV, the corrections are small. 


1 Introduction 


The search for the Higgs boson is one of the most important objectives of present and future 
colliders. A Higgs boson or some object like it is needed in order to give the and Z 
gauge bosons their observed masses and to cancel the divergences which arise when radiative 
corrections to electroweak observables are computed. We have few clues as to the expected 
mass of the Higgs boson, which a priori is a free parameter of the theory. Direct experimental 
searches for the Standard Model Higgs boson at LEP and LEP2 yield the limit, 

Mh > 89.8 GeV at 95% c.l. (1) 

LEP2 will eventually extend this limit to around > 109 GeV at 95% c.l. Also, from the 
analysis of the electroweak precision measurements an indirect upper bound can be found, [Q] 

Mh < 280 GeV at 95% c.l. ( 2 ) 

Above the LEP2 direct search limit and below the 2Z boson pair production threshold 
is termed the intermediate mass region and is the most difficult Higgs mass region to probe 
experimentally. In this mass range the associated production of a Higgs boson with either a 
gauge boson or a fermion - antifermion pair can be an important discovery channel. In this 
paper, we focus on the associated production of a Higgs boson with a heavy quark pair. 

The associated production of a Higgs boson with a top quark pair in collisions has 

a small rate, around 1 fb for y/s = 500 GeV and Mh ~ 100 GeV. However, the signature, 
e’''e“ —>■ tth —>■ W^W~bbbb, is distinctive. [Q] The experimental viability of this signature has 
not yet been carefully evaluated. 

Once the Higgs boson has been discovered, it will be important to measure its couplings to 
fermions and gauge bosons. These couplings are completely determined in the Standard Model 
with no adjustable parameters. The couplings to the gauge bosons can be measured through 
the associated production processes, > Zh, qq' —>• W~^h, and qq —> Zh, and through 

vector boson fusion, W~^W~ —> h and ZZ ^ h. The couplings of the Higgs boson to fermions 
are more difficult to measure, however. |Q] 

The process e’''e“ ^ tth provides a direct mechanism for measuring the tth Yukawa coupling. 
Since this coupling can be significantly different in a supersymmetric model from that in the 
Standard Model, the measurement would provide a mean of discriminating between models. 
The tth Yukawa coupling also enters into the rates for gg ^ h and h 77 , as these processes 
have large contributions from top quark loops. However, in these cases it is possible that there 
is unknown new physics which also enters into the rate and dilutes the interpretation of the 
signal as the measurement of the tth coupling. Ref. [Q estimates that for a Higgs boson with 
mass less than 100 GeV, both the tth and bbh couplings can be measured to an accuracy of 
roughly ±30 % at the LHG with 600 fb~^ of data. It is possible that a high energy lepton 
collider could obtain a higher precision on the measurement of the couplings. 
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The QCD corrections to the associated production of a Higgs boson with a heavy quark 
pair are the subject of this paper. They have been independently computed by Dittmaier 
et.al. @] We compute the QCD corrections to the process > tih, including only the 

photon exchange diagrams, which constitute the dominant contribution to the cross section. 
Our corrections are valid for all values of the Higgs boson mass and for any energy. Section ^ 
contains a review of the lowest order rate, including a discussion of the relative importance 
of the Z boson contribution (which we neglect) and the Higgs bremsstrahlung from the Z 
boson. The Higgs bremsstrahlung from the Z is always less than a few percent effect and so 
the e’''e“ —> tih process remains a candidate for measuring the coupling of the Higgs to the top 
quark. 

The 0{as) contributions are discussed in Sections and The contributions from the gluon 
bremsstrahlung diagrams and the separation of the calculation of the gluon emission into hard 
and soft pieces are presented in Section |^. The one-loop virtual contributions are discussed in 
Section ^ Our numerical results for a/s = 500 GeV and 1 TeV are presented in Section 

The corrections to e’'“e“ —tth have previously been estimated in an approximation valid 
at high energy and for Mh << Mt. [^. In the region where the approximation is valid, 
a/s > 1 TeV and < 100 GeV, the QGD corrections were estimated to be small. We end 
with a comparison of our results with the approximate calculation. 

2 e^e~ tth: Lowest Order 

e~ t e~ t 


+ ■■■ 

t 

Figure 1: Feynman diagrams contributing to the lowest order process, e+e“ — tth. 

The cross section for tih occurs through the Feynman diagrams of Fig. |l| and 

was hrst calculated in ^ (photon-exchange contribution only) and then completed in 0 (both 
photon and Z-exchange contributions). We write the cross section in the form 

da{e'^e~ —> tih^) 
dxh 
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Figure 2: Lowest order cross section for tth at ^/s = 500 GeV and y/s = 1 TeV. The 

curve labelled complete includes both 7 and Z exchange, along with bremsstrahlung from the 
Z boson. We take = 175 GeV. 

where (Tq = dvra^/Ss, a is the QED hne structure constant, W = 3 is the number of colors, 
Xh = 2Eh/^/s with Eh the Higgs boson energy, and Qi, V and Ai (i = e, t) denote the 
electromagnetic and weak couplings of the electron and of the top quark respectively, 

V- = ~ 4QiS^ ^ 211^ 

4swCiv ’ * 4swCw 

with = ±1/2 being the weak isospin of the left-handed fermions and = 1 — = 0.23. 

(The contribution from the photon alone can be trivially found by setting V = = 0). 

The coefficients Gi and G 2 describe the radiation of the Higgs boson off the top quark (both 
photon and Z boson exchange) and are given by. 


Gi = 




(/52 - xl) Xh 


4/3 (4M/ - Ml) (2Ml + s) Xh+ 


( 5 ) 


- xi) (I 6 M/ ± 2Ml - 2MlsXh + s^xl - 4Ml (sMl -2s- 2sxh)) log 

\Xh P, 
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where gt is the top quark Yukawa coupling, 
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( 7 ) 


with = 2Mh/y/s and = 1 — AM^ /s + Ml/s. 

The other four coefficients, G^,... ,Gq describe the emission of a Higgs boson from the Z-boson 
and can be written in the following form. 


Ga 
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G 5 
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where gz denotes the coupling of the Higgs boson to the Z boson. 


Gfi = 


'Xh + 

yXh-i3, 

^Xh + 

KXh -/3. 


( 8 ) 


gz = 



(9) 


V 

As already observed in Ref. 0], the most relevant contributions are those in which the Higgs 
boson is emitted from a top quark leg, i.e. those proportional to Gi and G 2 in Eq. (|^). The 
contribution from the Higgs boson coupling to the Z boson is always less than a few per cent 
at y/s = 500 GeV and 1 TeV and can safely be neglected. In Fig. |], we show the complete 
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cross section for e’''e“ —^ tth production and also the contribution from the photon exchange 
contribution only. We see that at both a/s = 500 GeV and 1 TeV, the cross section is very well 
approximated by the photon exchange only. In the remainder of the paper, we will consider only 
the photon exchange contribution, neglecting the Z boson exchange contribution everywhere. 


3 Real Gluon Emission 

The 0{as) inclusive cross section for —>■ tth receives contributions from real gluon emission 

from the hnal quark legs, 

e+(p) + e"(g) ^ t{pt) + t{p't) + KPh) + 9 {k) , (10) 

as shown in Fig. The cross section can be separated into hard and soft pieces by introducing 
an arbitrary separation Emin on fhe gluon momenta, such that for Eg < Emin-, we can use the 
eikonal approximation and calculate the cross section analytically, while for Eg > Emin-, we 
integrate over the phase space numerically. The hnal result is of course independent of the 
cut-off Emin (see Section p.2|) . 

3.1 Soft Gluon Radiation 

For soft gluons. Eg < Emin) we neglect the momenta of the radiated gluons everywhere but in 
the singular propagators. In the soft gluon approximation, we hnd the contribution from the 
radiated gluons, 

(—\ =-(—\ f _ 2pt ■ pt' \ 

\dxh)soft \dxh)) J\k\<Ernm2.u}{k)\{pt-ky {pf ■ ky {pt ■ k){pt'■ k)} ' 

where k is the gluon momentum, a;(/c) = | fc / +mg and we have introduced a small gluon 

mass rUg to regulate the infrared divergences occurring in the soft gluon emission. The depen¬ 
dence on the gluon mass will be cancelled by contributions from the virtual graphs which are 
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Figure 4: Contributions to the process e+e“ —tthg as a function of the cut-off on the soft gluon 
energy, Emini for hxed gluon and Higgs boson masses at a/s = 500 GeV. We take Mt = 175 GeV 
and as{Mt) = .11164. 


also evaluated with a non-zero gluon mass (see Section [4.3|) . The lowest order result can be 
found from Eq. (||). 

The integral over the soft gluon phase space has been performed in Refs. ^ and 0 leading 
to the analytic result for the soft contribution to e^e“ —>■ tthg, 
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where x is the solution to 
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Figure 5: Contributions from the process e^e“ —>■ tthg as a function of the cut-off on the 
soft gluon energy, -Emm, for fixed gluon and Higgs boson masses at a/s = 1 TeV. We take 
Mt = 175 GeV and as{Mt) = .11164 


subject to the constraint 


+ 1) - 2xpt ■ pt' = 0 , 


and 


- Pt' 
„o 
Pt' 


> 0 




- 1 ) 
2{xp^t-Pt') 


(13) 


(14) 


(15) 


3.2 Hard gluon radiation 

The hard gluon contribution is calculated from the diagrams of Fig. ^ for gluon momenta 
Eg > Emin- The integration over the final state phase space is done using a numerical Monte 
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Figure 6: Box diagrams contributing to the virtual corrections to e’''e ^ tth. 


Carlo. The hard gluon cross section is insensitive to the small value chosen for the gluon mass. 
In Figs. ^ and ^ we show the contributions of the soft and hard gluon radiation as a function 
of Emin- The sum of the soft and hard terms is clearly independent of the separation, Emin 
for Emin >> iT^g- The result retains a dependence on the gluon mass, rrig, however, which 
is cancelled when the infrared divergent pieces of the virtual contributions are included (see 
Section ||). 

We choose v = ~ 246 GeV, Mt = 175 GeV, and use the 1-loop value for 

as{Mt) = .11164 (which corresponds to Aqcd = 200 MeV) in all our numerical calculations. 


4 Virtual Corrections 

The amplitude for e’''e“ —>■ tth including virtual corrections to 0{as) can be written as 

( 1 / 

^ , (16) 
47r 

where Cp = 4/3. We include in also the wave function, Yukawa coupling, and mass 
counterterms. The corresponding contribution to the cross section at 0{as) is 

^virt _ 2Re (^f U*) = ^G^Re (^f M*) . (17) 



Figure 7: Vertex diagrams of type 1 contributing to the virtual corrections to e'''e 


tth. 




Figure 9: Feynman diagrams contributing to the self energy corrections of the internal heavy 
quark propagators for tth. 


The virtual corrections we consider are: 


• box diagrams (2 diagrams, see Fig. ||); 

• vertex corrections of type 1 (2 diagrams, see Fig. ^ ); 

• vertex corrections of type 2 (2 diagrams, see Fig. 

• self-energy corrections to the internal leg propagators (2 diagrams, see Fig. |^); 

• wave function and Yukawa coupling renormalization; 

• mass renormalization (2 counterterm insertions, one for each internal propagator). 


In general, virtual diagrams can have both ultraviolet (UV) and infrared (IR) singularities. 
In order to explicitly check the cancellation of both UV and IR divergences among different 
loop diagrams and counterterms, we have regularized them using dimensional regularization. 
We have then verihed the cancellation of the 1/euv and 1/eiR poles both analytically and 
numerically. We will discuss IR divergences in greater detail in Section 4.3. 
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Figure 10: Contributions to e’''e —tth from virtual diagrams at a/s = 500 GeV with a finite 
gluon mass. We take Mt = 175 GeV and as(M/ = .11164 

4.1 Wave Function and Mass Renormalization Counterterms 

The wave function renormalization contribution to the cross section can be expressed directly 
in terms of Uo, the cross section at the tree level for —>• tih, 

<^^2 = . (18) 

The tree level cross section must be computed to C>(e), hence we denote it by (Tq. We separate 
the infrared and ultraviolet divergences and write the wave function renormalization countertern 
af|, 


= 1 + 5Z2 


(19) 


^We thank the authors of Ref. [4] for pointing out an inconsistency in our original treatment of the wave- 
function renormalization. 
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The renormalization of the Yukawa coupling constant contributes to the 0(as) cross section, 


where we use the pole dehnition of the top quark mass, 

- ^^— 2 ' 


Lt _ ots 
Mt ~ 




( 20 ) 


( 21 ) 


Finally, there is the mass renormalization of the internal heavy quark propagators which is 
calculated using Eq. ( 0 )- 


4.2 Calculation of Loop Diagrams 

The analytical reduction of the one loop diagrams to a linear combination of one loop tensor 
and scalar integrals is performed using FORM, MAPLE and MAXIMA. The one loop integrals 
are then evaluated with the help of the numerical program FF [^]. For a hxed set of momenta. 


this program evaluates the hnite contribution to the loop integrals numerically. The integrals 
over the phase space are performed with a numerical Monte Carlo. In Fig. [l^, we show the 
finite contributions from the virtual diagrams. The virtual contribution can be written as 


(^virt (^box T CTyi T T (^self T count i (22) 

where abox, o'vi, crv 2 and ageif are shown in Figs. |^ - The counterterms are acount = o'z 2 +o'y + 
am, where am is the contribution of the internal propagator mass renormalization. Individually, 
the diagrams generate \og{^^/M^) contributions which are associated with 1/euv singularities 
and we take the renormalization scale /i = Mt. Some virtual corrections also generate IR 
divergences and we discuss these in Section |4.3| . The large cancellations between the various 
contributions is clear from Fig. 0 


4.3 Cancellation of Infrared Divergences 

As we discussed before, we have verified the cancellation of the 1/euv ultraviolet singularities 
both analytically and numerically. The cancellation of the l/em singularities has also been 
verihed analytically. Equivalently, the infrared divergences can be regulated by introducing 
a hnite gluon mass, as was done for the real gluon diagrams of Section This amounts to 
making the substitution 


F(e)(47r/i2)" ^ log(mp 
which is what we have used in our numerical calculations. 


(23) 
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Figure 11: Contributions to e^e —> tth which have infrared divergences. We take Mt = 
175 GeV and as{Mt) = .11164 

Infrared divergences arise from the box diagram, Fig. P, the Z 2 wavefunction renormaliza¬ 
tion, Eq. (p0[), and the soft gluon emission, Eq. ([T^). In Figs. and we show the sum of 
these contributions as a function of the gluon mass and see that the sum is independent of the 
gluon mass, conhrming the cancellation of the infrared divergences. 


5 Results 


5.1 Numerical Results 


When the virtual and real contributions are combined, the hnal result is hnite and independent 
of both Emin and rUg. In Figs. and we show the various contributions to the total cross 
section, ui is the complete 0{as) corrected rate. 


4” ^virt 4“ ^hard 4“ ^soft ■ 

The combination ayirt 4- Usoft is independent of the gluon mass, but retains a dependence on 
Emin which is cancelled by (Thard- At ^/s = 500 GeV, the corrections are large and positive. 
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Figure 12: Contributions to e^e —> tth which have infrared divergences. We take Mt = 
175 GeV and as{Mt) = .11164 

signihcantly increasing the rate. The corrections are smaller at ^/s = 1 TeV, with large cancel¬ 
lations between the hard and the virtual plus soft contributions. 

The size of the QCD corrections can be described by a iC factor, 

KM = - . ( 25 ) 

which is shown in Figs. and 0. Note that after the cancellation of the 1/euv divergences, 
the only p dependence is in If /i = ^/s, then K{Mh = 100 GeV) is reduced to 1.4 from 

the value K = 1.5 obtained with jj, = Mf for a/s = 500 GeV. The authors of Ref. choose 
/i = a/s. Our results for the K factor are in agreement with theirs. 

In Fig. 0 we show the shape of the differential cross section ^ at ^/s = 500 GeV and for 
Mh = 100 GeV. The cross section is peaked around Xh ~ -45. Including the 0{as) corrections 
has little effect on the shape of the distribution. 
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Figure 13: QCD corrections to e’''e —> tth at a/s = 500 GeV. cxo is the lowest order cross section 
and cTi is the complete 0{as) corrected rate. We take Mf = 175 GeV and aa{Mt) = .11164 

5.2 Comparison with Approximate Result 

In Ref. [Q, the 0{as) corrections to the process e^e“ —>■ tth are computed in a framework 
where the Higgs boson is treated as a parton which is radiated from the heavy top quark. The 
distribution of Higgs bosons in a top quark, ft-,h{xh) is computed to 0{as) and convoluted 
with the cross section for ^ tt, also computed to 0{as), 

^max 

cr(e+e" ^ tth)EHA = 2 / ^ dxhft^h{xh)a{e^e~ tt) . (26) 

} rrxnin 

This approximation, which we call the Effective Higgs Approximation (EHA), is expected to 
be valid for Mh « Mt « \/s. 

The impact of QGD corrections on the prediction of the cross section for e’''e“ ^ tih is 
described by the K-factor of Eq. (pS]), where the cross sections are evaluated using Eq. (p^) 
to the appropriate order in as- In Ref. 0, we found = a/s) — -94 at a/s = 1 TeV 
for Mh ~ 100 GeV. Gomparing with Fig. 0 we see good agreement between the EHA at 
a/s = 1 TeV and the exact calculation presented here. However, the hard gluon bremsstrahlung 
terms of Section cannot be adequately included within the context of the EHA. From Fig. ^ 
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Figure 14: QCD corrections to e'''e —> tth at a/s = 1 TeV. ao is the lowest order cross section 
and cTi is the complete 0{as) corrected rate. We take Mf = 175 GeV and aa{Mt) = .11164 

it is apparent that the hard gluon terms are not small. Hence the agreement between the EHA 
and the present calculation seems to derive from some more complicated cancellation between 
the hard gluon terms and those 0{Mf /s), 0{Ml/s) and 0{Ml/M^) contributions that are also 
sistematically neglected in the EHA. 

6 Conclusion 

We have computed the 0{as) corrected rate for tth. At a/s = 500 GeV, the corrections 

are large and positive, while at a/s = 1 TeV, the QGD corrections are small. Studies of the 
experimental viability of this process as a means of measuring the top quark -Higgs boson 
Yukawa coupling are needed in order to assess the usefulness of the tth process. 
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Figure 15: Ratio of the 0(as) corrected rate to the lowest order cross section for e’''e —> tth 
at a/s = 500 GeV. We take Mt = 175 GeV and as{Mt) = .11164 
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Figure 17: Differential cross section, da/dxh, for —>• ith for Mh = 100 GeV and = 
500 GeV. (Jo is the lowest order cross section and ai is the complete 0{as) corrected rate. We 
take Mt = 175 GeV and as{Mt) = .11164. 
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